The effect of ambipolar electric fields on the electron heating in capacitive RF plasmas

We investigate the electron heating dynamics in electropositive argon and helium capacitively coupled RF discharges driven at 13.56 MHz by Particle in Cell simulations and by an analytical model. The model allows to calculate the electric field outside the electrode sheaths, space and time resolved within the RF period. Electrons are found to be heated by strong ambipolar electric fields outside the sheath during the phase of sheath expansion in addition to classical sheath expansion heating. By tracing individual electrons we also show that ionization is primarily caused by electrons that collide with the expanding sheath edge multiple times during one phase of sheath expansion due to backscattering towards the sheath by collisions. A synergistic combination of these different heating events during one phase of sheath expansion is required to accelerate an electron to energies above the threshold for ionization. The ambipolar electric field outside the sheath is found to be time modulated due to a time modulation of the electron mean energy caused by the presence of sheath expansion heating only during one half of the RF period at a given electrode. This time modulation results in more electron heating than cooling inside the region of high electric field outside the sheath on time average. If an electric field reversal is present during sheath collapse, this time modulation and, thus, the asymmetry between the phases of sheath expansion and collapse will be enhanced. We propose that the ambipolar electron heating should be included in models describing electron heating in capacitive RF plasmas

Customized ion flux-energy distribution functions in capacitively coupled plasmas by voltage waveform tailoring

We propose a method to generate a single peak at a distinct energy in the ion flux-energy distribution function (IDF) at the electrode surfaces in capacitively coupled plasmas. The technique is based on the tailoring of the driving voltage waveform, i.e. adjusting the phases and amplitudes of the applied harmonics, to optimize the accumulation of ions created by charge exchange collisions and their subsequent acceleration by the sheath electric field. The position of the peak (i.e. the ion energy) and the flux of the ions within the peak of the IDF can be controlled in a wide domain by tuning the parameters of the applied RF voltage waveform, allowing optimization of various applications where surface reactions are induced at particular ion energies

OPerationalized ALignment: Assessing feasibility of value constellations exploiting innovative services

Tan, Y.H. [Promotor]Akkermans, J.M. [Copromotor

Effects of fast atoms and energy-dependent secondary electron emission yields in PIC/MCC simulations of capacitively coupled plasmas

In most PIC/MCC simulations of radio frequency capacitively coupled plasmas (CCPs) several simplifications are made: (i) fast neutrals are not traced, (ii) heavy particle induced excitation and ionization are neglected, (iii) secondary electron emission from boundary surfaces due to neutral particle impact is not taken into account, and (iv) the secondary electron emission coefficient is assumed to be constant, i.e. independent of the incident particle energy and the surface conditions. Here we question the validity of these simplifications under conditions typical for plasma processing applications. We study the effects of including fast neutrals and using realistic energy-dependent secondary electron emission coefficients for ions and fast neutrals in simulations of CCPs operated in argon at 13.56 MHz and at neutral gas pressures between 3 Pa and 100 Pa. We find a strong increase of the plasma density and the ion flux to the electrodes under most conditions, if these processes are included realistically in the simulation. The sheath widths are found to be significantly smaller and the simulation is found to diverge at high pressures for high voltage amplitudes in qualitative agreement with experimental findings. By switching individual processes on and off in the simulation we identify their individual effects on the ionization dynamics and plasma parameters. We conclude that fast neutrals and energy-dependent secondary electron emission coefficients must be included in simulations of CCPs in order to yield realistic results

Kinetic simulation of the sheath dynamics in the intermediate radio-frequency regime

The dynamics of temporally modulated plasma boundary sheaths is studied in the intermediate radio frequency regime where the applied radio frequency and the ion plasma frequency are comparable. Two kinetic simulation codes are employed and their results are compared. The first code is a realization of the well-known scheme, Particle-In-Cell with Monte Carlo collisions (PIC/MCC) and simulates the entire discharge, a planar radio frequency capacitively coupled plasma (RF-CCP) with an additional heating source. The second code is based on the recently published scheme Ensemble-in-Spacetime (EST); it resolves only the sheath and requires the time resolved voltage across and the ion flux into the sheath as input. Ion inertia causes a temporal asymmetry (hysteresis) of the sheath charge-voltage relation; also other ion transit time effects are found. The two codes are in good agreement, both with respect to the spatial and temporal dynamics of the sheath and with respect to the ion energy distributions at the electrodes. It is concluded that the EST scheme may serve as an efficient post-processor for fluid or global simulations and for measurements: It can rapidly and accurately calculate ion distribution functions even when no genuine kinetic information is available

Kinetic Interpretation of Resonance Phenomena in Low Pressure Capacitively Coupled Radio Frequency Plasmas

The kinetic origin of resonance phenomena in capacitively coupled radio frequency plasmas is discovered based on particle-based numerical simulations. The analysis of the spatio-temporal distributions of plasma parameters such as the densities of hot and cold electrons, as well as the conduction and displacement currents reveals the mechanism of the formation of multiple electron beams during sheath expansion. The interplay between highly energetic beam electrons and low energetic bulk electrons is identified as the physical origin of the excitation of harmonics in the current

Simulation benchmarks for low-pressure plasmas: capacitive discharges

Benchmarking is generally accepted as an important element in demonstrating the correctness of computer simulations. In the modern sense, a benchmark is a computer simulation result that has evidence of correctness, is accompanied by estimates of relevant errors, and which can thus be used as a basis for judging the accuracy and efficiency of other codes. In this paper, we present four benchmark cases related to capacitively coupled discharges. These benchmarks prescribe all relevant physical and numerical parameters. We have simulated the benchmark conditions using five independently developed particle-in-cell codes. We show that the results of these simulations are statistically indistinguishable, within bounds of uncertainty that we define. We therefore claim that the results of these simulations represent strong benchmarks, that can be used as a basis for evaluating the accuracy of other codes. These other codes could include other approaches than particle-in-cell simulations, where benchmarking could examine not just implementation accuracy and efficiency, but also the fidelity of different physical models, such as moment or hybrid models. We discuss an example of this kind in an appendix. Of course, the methodology that we have developed can also be readily extended to a suite of benchmarks with coverage of a wider range of physical and chemical phenomena

The effect of the driving frequency on the confinement of beam electrons and plasma density in low pressure capacitive discharges

The effect of changing the driving frequency on the plasma density and the electron dynamics in a capacitive radio-frequency argon plasma operated at low pressures of a few Pa is investigated by Particle in Cell/Monte Carlo Collisions simulations and analytical modeling. In contrast to previous assumptions the plasma density does not follow a quadratic dependence on the driving frequency in this non-local collisionless regime. Instead, a step-like increase at a distinct driving frequency is observed. Based on the analytical power balance model, in combination with a detailed analysis of the electron kinetics, the density jump is found to be caused by an electron heating mode transition from the classical $\alpha$-mode into a low density resonant heating mode characterized by the generation of two energetic electron beams at each electrode per sheath expansion phase. These electron beams propagate through the bulk without collisions and interact with the opposing sheath. In the low density mode, the second beam is found to hit the opposing sheath during its collapse. Consequently, a high number of energetic electrons is lost at the electrodes resulting in a poor confinement of beam electrons in contrast to the classical $\alpha$-mode observed at higher driving frequencies. Based on the analytical model this modulated confinement quality and the related modulation of the energy lost per electron lost at the electrodes is demonstrated to cause the step-like change of the plasma density. The effects of a variation of the electrode gap, the neutral gas pressure, the electron sticking and secondary electron emission coefficients of the electrodes on this step-like increase of the plasma density are analyzed based on the simulation results

Tailored voltage waveform capacitively coupled plasmas in electronegative gases : frequency dependence of asymmetry effects

Capacitively coupled radio frequency plasmas operated in an electronegative gas (CF4) and driven by voltage waveforms composed of four consecutive harmonics are investigated for different fundamental driving frequencies using PIC/MCC simulations and an analytical model. As has been observed previously for electropositive gases, the application of peak-shaped waveforms (that are characterized by a strong amplitude asymmetry) results in the development of a DC self-bias due to the electrical asymmetry effect (EAE), which increases the energy of ions arriving at the powered electrode. In contrast to the electropositive case (Korolov et al 2012 J. Phys. D: Appl. Phys. 45 465202) the absolute value of the DC self-bias is found to increase as the fundamental frequency is reduced in this electronegative discharge, providing an increased range over which the DC self-bias can be controlled. The analytical model reveals that this increased DC self-bias is caused by changes in the spatial profile and the mean value of the net charge density in the grounded electrode sheath. The spatio-temporally resolved simulation data show that as the frequency is reduced the grounded electrode sheath region becomes electronegative. The presence of negative ions in this sheath leads to very different dynamics of the power absorption of electrons, which in turn enhances the local electronegativity and plasma density via ionization and attachment processes. The ion flux to the grounded electrode (where the ion energy is lowest) can be up to twice that to the powered electrode. At the same time, while the mean ion energies at both electrodes are quite different, their ratio remains approximately constant for all base frequencies studied here